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Dilepton and photon production in heavy-ion collisions at SPS energies are studied in 
the relativistic transport model that incorporates self-consistently the change of hadron 
masses in dense matter. It is found that the dilepton spectra in proton-nucleus reactions 
can be well described by the conventional mechanism of Dalitz decay and direct vector 
meson decay. However, to provide a quantitative explanation of the observed dilepton 
spectra in central heavy-ion collisions requires contributions other than these direct de¬ 
cays and also various medium effects. Introducing a decrease of vector meson masses 
in hot dense medium, we find that these heavy-ion data can be satisfactorily explained. 
Furthermore, the single photon spectra in our calculations with either free or in-medium 
meson masses do not exceed the upper bound deduced from the experiments by the WA80 
Collaboration. 


1. INTRODUCTION 


The measurement of electromagnetic observables such as the photon and dilepton spec¬ 
tra constitutes a major part of CERN-SPS heavy-ion programs. It will soon become one 
of major efforts at the BNL-RHIC heavy-ion programs. The primary reason for this is 
that photons and dileptons do not suffer strong final-state interactions as hadrons do. 
They can thus be considered as ‘penetrating probes’ of the initial hot and dense stages 
involved in high-energy heavy-ion collisions. 

Recent observation of the enhancement of low-mass dileptons in central heavy-ion col¬ 
lisions at CERN-SPS energies |Tj.j2|] has generated a great deal of interest in the heavy- 
ion community. Different dynamical models, such as the hydrodynamical and transport 
models, have been used to investigate this phenomenon |3| [T2||. Calculations based on 
‘conventional sources’ such as the Dalitz decay and direct vector meson decay that ac¬ 
count for the dilepton spectra in proton-induced reactions failed to explain the observed 
enhancement in heavy-ion collisions. Various medium effects, such as the dropping vector 
meson masses [H-El-lllIIl] as first proposed by Brown and Rho [|l3|], the modification of 
rho meson spectral function ||. an d the enhanced production of r] and/or rj' |I4|, have 
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been proposed to explain this enhancement. The last possibility has been ruled out by 
the measured g spectra in the same reactions (T^j. The effects due to medium modifi¬ 


cations of rho meson spectral function on low-mass dileptons have been examined using 
simplified fireball models [;6j; questions remain concerning the exact space-time evolution 
and baryon chemical compositions which play an important role in the determination of 
the rho meson spectral function. In this contribution we will concentrate on the idea of 
dropping vector meson masses, which is the focus of Section 3. 

Another piece of experimental data from CERN-SPS that has been discussed extensively 
is the single photon spectra from the WA80 collaboration |T6|. In various hydrodynamics 
calculations the absence of significant thermal photon production has been used as an 
indication of quark gluon plasma formation that lowers the initial temperature |L7|. How¬ 
ever, detailed transport model calculations find this conclusion to be premature. This 
will be discussed in Section 4. A brief summary is given in Section 5. 


2. THE RELATIVISTIC TRANSPORT MODEL 


In studying medium effects in heavy-ion collisions, the relativistic transport model |L8] 
based on the Walecka-type model |T9[ has been quite useful, as it provides a thermody¬ 
namically consistent description of the medium effects through the scalar and vector fields. 
In heavy-ion collisions at CERN-SPS energies, many hadrons are produced in the initial 
nucleon-nucleon interactions. This is usually modeled by the fragmentation of strings, 
which are the chromoelectric flux-tubes excited from the interacting quarks. One suc¬ 
cessful model for taking into account this nonequilibrium dynamics is the RQMD model 
20[]. To extend the relativistic transport model to heavy-ion collisions at these energies, 


we have used as initial conditions the hadron abundance and distributions obtained from 
the string fragmentation in RQMD. Further interactions and decays of these hadrons are 
then taken into account as in usual relativistic transport model. 

To study the effects of dropping vector meson masses 


13,21] on the dilepton spectrum 


in heavy-ion collisions, we have extended the Walecka model from the coupling of nucleons 
to scalar and vector fields to the coupling of light quarks to these fields, using the ideas 
of the meson-quark coupling model |22|. For a system of nucleons, pseudoscalar mesons, 


vector mesons, and axial-vector mesons at temperature T and baryon density ps, the 
scalar field (a) is determined self-consistently from 
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where we have used the constituent quark model relations for the nucleon and vector 
meson masses [|2zj, i.e., m* N = — g a (cr ), m* ~ m 

of the ij meson in free space which leads to m* ~ 
rule relation between the rho-meson and a\ meson masses, i.e, m 


- (2/3)g a (a ), the quark structure 
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We recently found that the use of a refined model, the effective chiral Lagrangian of |23 
led to essentially the same results for dilepton spectra p4 



Figure 1. Dilepton invariant mass spectra in p+Be (left window) and p+Au (right win¬ 
dow) collisions at 450 GeV. 


3. DILEPTON PRODUCTION 


The main contributions to dileptons with mass below 1.2 GeV are the Dalitz decay of 
7r°, r) and u, the direct leptonic decay of p°, lv and (j), the pion-pion annihilation which 
proceeds through the p° meson, and the kaon-antikaon annihilation that proceeds through 
the (j) meson. The differential widths for the Dalitz decay of 7T°, p, and u are related to 
their radiative decay widths via the vector dominance model, which are taken from Ref. 


25 


The decay of a vector meson into the dilepton is determined by the width, 


m 
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The coefficient Ci+i- in the dielectron channel is 8.814x 1CT 6 , 0.767x 10~ 6 , and 1.344x 10 -6 
for p, u, and 0, respectively, and is determined from the measured width. For the dimuon 
channel, these values are slightly larger. 
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In our model, dileptons are emitted continuously during the time evolution of the 
colliding system. The way the dilepton yield is calculated can be illustrated by the decay 
of a rho meson. Denoting, at time t, the differential multiplicity of neutral rho mesons by 
dN f) o (t) / dM , then the differential dilepton production probability is given by 


dN l+l - 

dM 


df dN p o(t) 
3 dM 


r p o. 


*i+i-(M)dt + 


dN p o ( t f ) TpQ(M) 
dM Tp(M) 


(3) 


where tf is the freeze-out time, which is found to be about 20 frn/c. The first term 
corresponds to dilepton emission before freeze out while the second term is from the 
decay of rho mesons still present at freeze out. 



M (MeV) M (MeV) 


Figure 2. Dilepton invariant mass spectra in p+W (left window) and S+W (right window) 
collisions at 200 AGeV. In the right window the solid and dashed curves are obtained with 
in-medium and free meson masses, respectively. 


The results for dilepton spectra from p+Be and p+Au collisions at 450 GeV are shown 
in Fig. 1, together with data from the CERES [|lj. It is seen that the data can be well 
reproduced by Dalitz decay of 7r°, rj and u mesons, and direct leptonic decay of p°, u and </> 
mesons. These results are thus similar to those found in Ref. [|4| using the Hadron-String 
Dynamics and those constructed by the CERES collaboration from known and expected 
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sources of clileptons [1]. Similar conclusions can be drawn concerning dimuon spectra in 
p+W collisions at 200 GeV, which are shown in the left window of Fig. 2. 

Our results for dilepton spectra in central S+Au collisions are shown in the left window 
of Fig. 3. The dashed curve is obtained with free meson masses. Although pion-pion 
annihilation is important for dileptons with invariant mass from 0.3 to 0.65 GeV, it still 
does not give enough number of clileptons in this mass region. Furthermore, for masses 
around m P}U there are more clileptons predicted by the theoretical calculations than shown 
in the experimental data. These are very similar to the results of Cassing et al. [fjj based on 
the Hadron-String Dynamics model and Srivastava et al. [j5| based on the hydrodynamical 
model. The results obtained with in-medium meson masses are shown by the solid curve. 
Compared with the results obtained with free meson masses, there is about a factor of 
2-3 enhancement of the dilepton yield in the mass region from 0.2 to 0.6 GeV, which thus 
leads to a good agreement between the theoretical results and the CERES data. Similar 
conclusions that dropping vector meson masses can explain the CERES dilepton data 
have been obtained in Refs. [HlHLin . 
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Figure 3. Dilepton invariant mass spectra in S+Au collision at 200 AGeV (left window) 
and Pb+Au collisions at 160 AGeV (right window). The solid and dashed curves are 
obtained with in-medium and free meson masses, respectively. 


In the right window of Fig. 3 we compare our prediction for central Au+Pb collisions at 
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160 AGeV with the preliminary data from the CERES collaboration. The normalization 
factor dN c hjdr / here is the average charge particle pseudo-rapidity density in the pseudo¬ 
rapidity range of 2 to 3, and is about 440 in this collision. In the results with free 
meson masses, shown by the dashed curve, there is a strong peak around m PtU , which 
is dominated by p° meson decay as a result of an enhanced contribution from pion-pion 
annihilation in Pb+Au collisions than in S+Au and proton-nucleus collisions. In the case 
of in-medium meson masses, shown by the solid curve, the p meson peak shifts to a lower 
mass, and the peak around m p ^ becomes a shoulder arising mainly from u meson decay. 
At the same time we see an enhancement of low-mass dileptons in the region of 0.25-0.6 
GeV as in S+Au collisions. The agreement with the data is thus significantly improved 
when dropping vector meson masses are used. 

The same model has been used to calculate the dimuon spectra from central S+W 
collisions. The results obtained with free meson masses are shown in the right window of 
Fig. 2 by the dashed curve, and are below the HELIOS-3 data in the mass region from 0.35 
to 0.6 GeV, and slightly above the data around m P:LJ as in the CERES case. However, the 
discrepancy between the theory and the data is somewhat smaller in this case due to the 
smaller charged-particle multiplicity at a larger rapidity than in the CERES experiment. 
Our results obtained with in-medium meson masses are shown in the right window of 
Fig. 2 by the solid curve, and are in good agreement with the data. The importance 
of dropping rho meson mass in explaining the HELIOS-3 data has also been found by 
Cassing et al. [Q. 


4. PHOTON PRODUCTION 


Single photon spectra in heavy-ion collisions at CERN-SPS energies have been mea¬ 
sured by the WA80 collaboration [|T6| . So far, only the upper bound has been determined 


for the so-called ‘thermal’ photon spectra, which accounts for about 5% of the total ob¬ 
served single photon yield, which is dominated by neutral pion decay. In hydrodynamical 
calculations 0 , the absence of significant thermal photons has been interpreted as an 
evidence for the formation of a quark gluon plasma. Without phase transition, the initial 
temperature of the hadronic gas found in Ref. JT?j is about 400 MeV. This leads to a 
large number of thermal photons from hadronic interaction which is not observed experi¬ 
mentally. Including the phase transition, the initial temperature is lowered to about 200 
MeV B because of increased degrees of freedom. With a lower initial temperature, the 
thermal photon yield is thus reduced and is found to better agree with the WA80 data. 
However, a more recent analysis |26| has shown that if one includes all hadron resonances 
with masses below 2.5 GeV the initial temperature can also be lowered to about 200 MeV, 
and the WA80 photon data can then be explained without invoking the formation of a 
quark gluon plasma. This is basically also the conclusion of the detailed transport model 
analysis presented here. 

For the ‘thermal’ photon spectra we include the decay of p, u;, r/, and (i\ mesons, 
as well as two-body processes such as 7T7 t —> py and up —> 7ry. The decay width for 
p —> 7T7T7 is taken from the model of Ref. [^7]] , which describes well the measured width 
for p° —> 7T7T7. Both the uj and a\ radiative decay widths are proportional to | | 3 ||28[, with 

the coefficients determined from the measured width. For the two-body cross sections, 
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we use the results of [^9|, which do not include the contribution from an intermediate a\ 


meson. The latter has already been included in our model as a two-step process. 

In Fig. 4, our results for ‘thermal’ photon spectra are compared with the upper bound 
from the WA80 collaboration. It is seen that the u meson radiative decay plays an im¬ 
portant role for photons with transverse momenta above 300 MeV. In both the free and 
in-medium meson mass scenarios, our thermal photon yields do not exceed the experimen¬ 
tal upper bound. Although there are more rho and a\ mesons in the case of in-medium 
meson masses, their contributions to single photon yield do not increase significantly, 
simply because their radiative decay widths decrease as a result of reduced phase space. 



Figure 4. ‘Thermal’ photon spectra in central S+Au collisions at 200 AGeV using free 
(left window) and in-medium (right window) meson masses. 


5. SUMMARY AND OUTLOOK 

In summary, we have studied dilepton production from both proton-nucleus and nucleus- 
nucleus collisions using the relativistic transport model with initial conditions determined 
by string fragmentation from the initial stage of the RQMD model. It is found that the 
dilepton spectra in proton-nucleus reactions measured by the CERES and the HELIOS-3 
collaboration can be well understood in terms of conventional mechanisms of Dalitz decay 













and direct vector meson decay. For dilepton spectra in central heavy-ion collisions, these 
conventional mechanisms, however, fail to explain the data, especially in the low-mass 
region from about 0.25 to about 0.6 GeV in CERES experiments, and from 0.35 to 0.65 
GeV in HELIOS-3 experiments. Including the contribution from pion-pion annihilation, 
which is important in the mass region from 2m n to m PttJ , removes some of the discrepancy. 
But the theoretical prediction is still substantially below the data in the low mass region 
and somewhat above the data around m p>u . The theoretical results are brought into 
good agreement with the data when reduced in-medium vector meson masses are taken 
into account. We have also calculated, within the same dynamical model, the thermal 
photon spectra in central S+Au collisions. In both free meson masses and in-medium 
meson masses scenarios, our results do not exceed the upper bound deduced from the 
experiments by the WA80 collaboration. 
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